Counting statistics of single-electron transport in a quantum dot 
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We have measured the full counting statistics (FCS) of current fluctuations in a semiconductor 
quantum dot (QD) by real-time detection of single electron tunneUng with a quantum point con- 
tact (QPC). This method gives direct access to the distribution function of current fluctuations. 
Suppression of the second moment (related to the shot noise) and the third moment (related to 
the asymmetry of the distribution) in a tunable semiconductor QD is demonstrated experimentally. 
With this method we demonstrate the ability to measure very low current and noise levels. 

PACS numbers: 72.70.+m, 73.23.Hk, 73.63.Kv 



Current fluctuations in conductors have been exten- 
sively studied because they provide additional informa- 
tion compared to the average current, in particular for 
interacting systems Q . Shot noise measurements demon- 
strated the charge of quasiparticles in the fractional 
quantum Hall effect 2] and in superconductors 3] . How- 
ever, to perform such measurements for semiconductor 
quantum dots (QD) using conventional noise measure- 
ments techniques is very challenging. This is because of 
the very low currents and the corresponding low noise lev- 
els in these systems. Earlier experiments demonstrated 
the measurement of shot noise in non-tunable QDs 4, 5], 
but to our knowledge, no experiments have been reported 
in the literature in which the tunnel barriers, and thereby 
the coupling symmetry, could be controlled 6]. 

An alternative way to investigate current fluctuations, 
introduced by Levitov et al.^ is known as full counting 
statistics (FCS) This method relies on the evalua- 
tion of the probability distribution function of the num- 
ber of electrons transferred through a conductor within 
a given time period. In addition to the current and the 
shot noise, which are the first and second moments of this 
distribution, this method gives access to higher order mo- 
ments. Of particular interest is the third moment (skew- 
ness), which is due to breaking the time reversal sym- 
metry at finite current. Experimentally, few attempts 
to measure the third moment have been made in tunnel 
junctions 0. 

The most intuitive method for measuring the FCS of 
electron transport is to count electrons passing one by 
one through the conductor. Real-time detection of single 
electron transport has been experimentally investigated 
only very recently It is a challenging task since 

it requires a very sensitive, low-noise and non-invasive 
electrometer, as well as a high-bandwidth circuit. Several 
devices, such as the single electron transistor | 9L lin l and 
the quantum point contact (QPC) 12, llElOEl, 
have been demonstrated to have high enough sensitivity 
to detect single electrons in a QD. But, up to now, none of 



these experiments were able to extract the full counting 
statistics of electron transport. 

Here we report on the real-time detection of single elec- 
tron tunneling through a QD using a QPC as a charge 
detector. With this method, we can directly measure 
the distribution function of current fluctuations in the 
QD by counting electrons. To our knowledge, this is the 
first measurement of the full counting statistics for elec- 
trons in a solid state device. In addition, we can tune 
the coupling of the QD with both leads and measure the 
respective tunneling rates. We show experimentally the 
suppression of the second and third moments of the cur- 
rent fluctuations when the QD is symmetrically coupled 
to the leads. 

FigureCIa) shows the structure, fabricated on a GaAs- 
GaAlAs heterostructure containing a two-dimensional 
electron gas 34 nm below the surface (density 4.5 x 10^^ 
m~^, mobility 25 m^(Vs)~^). An atomic force micro- 
scope (AFM) was used to oxidize locally the surface, 
thereby defining depleted regions below the oxide lines 
The measurements were performed in a ^He/^He 
dilution refrigerator with an electron temperature of 
about 350 mK, as determined from the width of ther- 
mally broadened Coulomb blockade resonances |^]. The 
charging energy of the QD is 2.1 meV and the mean 
level spacing is 200 — 300 /ieV. The conductance of the 
QPC, Gqpc, was tuned close to 0.25 x e'^/h. We apply 
a dc bias voltage between source and drain of the QPC, 
Vqpc = 500 /iV, and measure the current through the 
QPC, Iqpci which depends on the number of electrons 
TV in the QD. 

In order to measure the current with a charge detec- 
tor, one has to avoid that electrons travel back and forth 
between the dot and one lead or to the other lead due 
to thermal fluctuations [Fig. Q^b)]. This is achieved by 
applying a large bias voltage between source and drain, 
i.e. I ± eV/2 — Ed\ > kpT, where Ed is the electro- 
chemical potential of the dot and V is the symmetrically 
applied bias, see Fig. ^a, c). An example of a time 
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FIG. 1: (a) AFM micrograph of the sample consisting of a 
QD connected to two contacts S and and a nearby QPC. 
Gl, G2 and P are lateral gates allowing the tuning of the tun- 
nel coupling to the source S, the coupling to the drain D, and 
the conductance of the QPC. Gl and G2 are also used to tune 
the number of electrons in the QD. A symmetric bias voltage 
V is applied between the source and the drain on the QD. 
(b-c) Scheme of the quantum dot in the case of equilibrium 
charge fluctuations (b), and non-equiUbrium charge fluctua- 
tions (c). (d) Time trace of the current measured through 
the QPC corresponding to fluctuations of the charge of the 
dot between N and A/" + 1 electrons. The arrows indicate 
transitions where an electron is entering the QD from the 
source lead, (e) Probability density of the times nn and Tout 
(see text) obtained from time traces similar to the one in (d). 
The lines correspond to the expected exponential dependence 
(see the text), where the tunneling rates are calculated from 

'^/'^S{D) — ^/^in(out) — {^in{out))' 



trace of the QPC current in this configuration is shown 
in Fig. ^d). The number of electrons in the QD fluctu- 
ates between N and N^l. Since this trace corresponds to 
the non-equilibrium regime, we can attribute each tran- 
sition ^ A/" + 1 to an electron entering the QD from 
the source contact, and each transition A'^ 1 ^ A" to an 
electron leaving the QD to the drain contact. The charge 
fluctuations in the QD correspond to a non-equilibrium 
process, and are directly related to the current through 
the dot 10]. Due to Coulomb blockade, only one elec- 
tron at a time can enter the QD, which allows to count 
electrons traveling through the system. 

The first application of electron counting in the non- 
equilibrium regime concerns the determination of the in- 
dividual tunneling rates from the source to the QD, T^, 
and from the QD to the drain. Yd- Previous experiments 
determining the individual tunneling rates involved more 
than two leads connected to the QD ^19]. In the trace of 
Fig.njd), the time Tin corresponds to the time it takes for 
an electron to enter the QD from the source contact, and 
Tout to the time it takes for the electron to leave the QD to 
the drain contact. For independent tunneling events, the 
tunneling rates can be calculated from the average of Tin 
and Tout on a long time trace Q, ^/^s{d) = ^/^in{out) = 
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FIG. 2: Statistical distribution of the number n of electrons 
entering the QD during a given time to . The two panels corre- 
spond to two different values of the tunneling rates, obtained 
for different values of the gate voltage Vgi • The time to is cho- 
sen in order to have the same mean value of number of events, 
(n) ^ 3, for both graphs. We have checked that this choice 
does not affect the results. The line shows the theoretical 
distribution calculated from Eqs. (QJ and (0). The tunneling 
rates are determined experimentally by the method described 
in Fig.QJe), and no fitting parameters have been used for the 
theoretical curves. 



{rin{out))' To check that the tunneling events are indeed 
independent, we have compared the probability densi- 
ties and Prout with the expected exponential behav- 
ior p(r^n(ont)) = ^siD)^^p{-^siD)riniout))' FigurelUe) 
shows good agreement with our data. It is interesting 
to note that, in the case shown in Fig. de), the QD is 
almost symmetrically coupled. We demonstrate here a 
very sensitive method to determine the symmetry of the 
coupling alternative to Ref. 20. 

From traces similar to the one shown in Fig. ^d), we 
can directly determine the statistical properties of se- 
quential electron transport through the QD. We count 
the number n of electrons entering the QD from the 
source contact during a time period to, i.e. the num- 
ber of down-steps in Fig. ^d) (see arrows). We obtain 
the distribution function of n by repeating this counting 
procedure on m = T/to independent traces, T = 0.5 s 
being the total length of the time trace. The resulting 
distribution functions are shown for two different values 
of the tunneling rates in Figs.El^a) andEJb). 

The FCS theory allows to determine the distribution 
function of the number n of electrons traveling through 
a conductor 0: 
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where S{x) is the generating function, which has been 
calculated for a single level QD for large bias voltage 
\±eV/2-Ed\:>kBT ^: 
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Here F^ and Td are the effective tunneling rates, which 
take into account any possible spin degeneracy of the 
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FIG. 3: (a) Average number of electrons entering the QD, 
/i, measured as a function of the gate voltage Vgi and the 
bias voltage V. Far from the edges of the Coulomb blockade 
region, i.e. for | zb eV/2 — Ed\ ^ ksT, the fluctuations of n 
are directly related to current fluctuations. The dashed line 
correspond to the cross-section shown in panel (b) . (b) Three 
first moments of the fluctuations of n as a function of the 
bias voltage V and at a given gate voltage Vgi = —44 mV. 
The ground state (GS) as well as two excited states (ES) are 
clearly visible. The moments are scaled so that /i corresponds 
to the number of electrons entering the QD per second. In 
the gray region, the condition | ± eV/2 — Ed\ ^ ksT is not 
valid, and the number of electrons entering the QD cannot 
be taken as the current flowing through the QD. The width 
of this region is 9 x ksT/e ~ 300 /iV, determined from the 
width for which the Fermi distribution is between 0.01 and 
0.99. (c) Normalized second and third moments as a function 
of the bias voltage V and at a given gate voltage Vgi = —44 
mV. 



levels in the QD, and correspond to the tunneling rates 
we determine experimentally. We have calculated the 
theoretical distribution functions for the tunneling rates 
measured in the cases of Figs.E^a) andEl^b) [solid lines]. 
The agreement with the experimental distribution is very 
good, in particular, given that no fitting parameters were 
used. Both graphs show a clear qualitative difference: 
Figure I2Ib) shows a broader and more asymmetric dis- 
tribution than Fig. ISfa). We will see later that this dif- 
ference comes from the different asymmetries of the tun- 
neling rates. 

In order to perform a more quantitative analysis, we 
calculate the three first central moments given by /i = 
(n), and jii = ((n — (n))*) for i =2,3, where (...) rep- 
resents the mean over T/to periods of length to- The 
first moment (mean) gives access to the mean current, 
/ = e/i/to, and the second central moment (variance) to 
the shot noise. Si = 2e^/i2/^o (for much larger than 
the correlation time). We are also interested in the third 
central moment, /is, which gives the asymmetry of the 
distribution function around its maximum (skewness). 
An important difference to previous measurements of the 



third cumulant is that our method can be used to extract 
any higher order cumulants. For the data presented here, 
the accuracy of the higher cumulants is limited by the 
short length of the time traces. 

We first focus on the mean /i of the distribution. By 
measuring /i as a function of the voltage applied on gate 
Gl and the bias voltage V, we can construct the so-called 
Coulomb diamonds (see Fig. Ol^a)). The Coulomb dia- 
monds describe the charge stability of the QD, normally 
measured in standard transport experiments 6]. Here, 
we present a novel way of measuring Coulomb block- 
ade diamonds by time-resolved detection of the electrons 
using a non-invasive charge detector. We observe clear 
Coulomb blockade regions as well as regions of finite cur- 
rent. As we increase the bias voltage, we see steps in 
the current. The first step in Fig. IH^b) (see left arrow) 
corresponds to the alignment of the chemical potential 
of the source contact with the ground state in the QD, 
and the following steps with excited states in the QD. 
From the resolution of the Coulomb diamonds, we see 
that the sample is stable enough such that background 
charge fluctuations do not play a significant role 22]. 

In addition to the mean, we have calculated the sec- 
ond and third central moments of the electron counting 
statistics. These two moments are shown in Fig.|3Jb) for 
Vgi = —44 mV as a function of the bias voltage. The sec- 
ond moment (blue dotted line) reproduces the steps seen 
in the current. These two moments can be represented 
by their reduced quantities /i2/M (known as the Fano 
factor) and /isZ/i, as shown in Fig. [Sfc). Both normal- 
ized moments are almost independent of the bias voltage, 
and correspond to a reduction compared to the values 
^2/^ = ^3/^=1 expected for classical fluctuations with 
Poissonian counting statistics. Super-poissonian noise 
23] is not expected in our configuration. 

In a QD, one expects a reduction of the moments due 
to the fact that when one electron occupies the QD, a 
second electron cannot enter. This leads to correlations 
in the current fluctuations, and to a reduction of the 
noise. The reduction is maximal when the tunnel barriers 
are symmetric. For an asymmetrically coupled QD, the 
transport is governed by the slow barrier and the noise 
recovers the value for a single tunneling barrier. The 
normalized moments for a single level QD at large bias 
voltage can be expressed as a function of the asymmetry 
of the tunneling rates, a = {V s - T d) / (T s d) |21|]: 

^ = -(l + a^) and ^ = 1 (l + 3a^) . (3) 

The second central moment recovers earlier calculations 
of the Fano factor in a QD . We see in these equations 
that both moments are reduced for a symmetrically cou- 
pled QD (i.e. a = 0), and tend to the Poissonian values 
for an asymmetrically coupled QD (i.e. a = ±1). 

Reduction of the second moment (shot noise) due to 
Coulomb blockade has already been reported in the case 
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FIG. 4: (a) Second and (b) third normalized central moments 
of the fluctuations of n as a function of the asymmetry of the 
tunneling rates, a — (F^ — Vu)l{rs ^^d)- To increase the 
resolution, each point at a given asymmetry is obtained by 
averaging over about 50 points at a given voltage Vbi and 
in a window of bias voltage 1.5 < V < 3 mV. Error bars 
correspond to the standard error of this averaging process, 
and are of the size of the points if not shown. The dashed 
lines are the theoretical predictions given by Eqs. ©• No 
fitting parameters have been used, since the tunneling rates 
are fully determined experimentally (see Fig. ^e) and text). 
Inset of (b): Variation of the asymmetry of the tunneling 
rates, a, as a function of Vgi- 



of asymmetrically coupled QDs 0, Hf. In these exper- 
iments, reduction of the shot noise occurs due to bias 
voltage dependent effective tunneling rates [2^ . Here we 
report the reduction of the second, as well as the third 
moment for a fully controllable QD. In particular, we 
are able to continuously change the tunneling rates: by 
changing the gate voltage Vgi, we change the chemical 
potential in the QD, and also the asymmetry of the cou- 
pling by changing the opening of the source lead. The 
tunneling rates can be directly measured as described in 
Fig.^e), and the inset of Fig. ^b) shows the variation of 
asymmetry with gate voltage in the region of interest. In 
Fig. Efa) andElJb), we show the normalized second and 
third central moments as a function of the asymmetry a. 
The experimental data follow the theoretical predictions 
given by Eqs. very well. We note in particular that 
no fitting parameters have been used since the tunneling 
rates are determined experimentally. 

Our ability to measure the counting statistics of elec- 
tron transport relies on the high sensitivity of the QPC as 
a charge detector. The counting process that we demon- 
strate in this paper was not possible in previous exper- 
iments with the accuracy required for performing a sta- 
tistical analysis 10]. Given the bandwidth of our ex- 
perimental setup. A/ = 30 kHz, the method allows to 
measure currents up to 5 fA, and we can measure cur- 
rents as low as a few electrons per second, i.e., less than 
1 aA. The low-current limitation is mainly given by the 
length of the time trace and the stability of the QD, and 
is well below what can be measured with conventional 
current meters. In addition, as we directly count elec- 
trons one by one, this measurement is not sensitive to 
the noise and drifts of the experimental setup. It is also 
an very sensitive way of measuring low current noise lev- 



els. Conventional measurement techniques are usually 
limited by the current noise of the amplifiers (typically 
10-29 AVHz) iai: here we demonstrate a measure- 
ment of the noise power with a sensitivity better than 
10-3^ AVHz. 

In conclusion, we have measured current fluctuations 
in a semiconductor QD, using a QPC to detect single 
electron traveling through the QD. We show experimen- 
tally the reduction of the second and third moment of the 
distribution when the QD is symmetrically coupled to the 
leads. This ability to measure current fluctuations in a 
QD, as well as the very low noise level we demonstrate 
here, open new possibilities towards measuring electronic 
entanglement in quantum dot systems j25L [2a|. 
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